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Abstract 



The electrostatic potential generated by a point charge at rest in a simple static, 
spherically symmetric wormhole is given in the form of series of multipoles and in closed 
form. The general potential which is physically acceptable depends on a parameter due to 
the fact that the monopole solution is arbitrary. When the wormhole has ^-symmetry, 
the potential is completly determined. The calculation of the electrostatic self-energy and 
of the self-force is performed in all cases considered. 



^ . 1 Introduction 

O 

£0 • The determination in the form of series of multipoles of the electrostatic potential generated 

by a point charge at rest in some static, spherically symmetric wormholes has been a subject 
C^- 1 of investigation recently. The main purpose of all these works PJE] is to determine the 

electrostatic self-force acting on the point charge within this multipole formalism. We note 
that the self-force for the scalar field has been also considered [HE]. 

The determination in closed form of the electrostatic potential is interesting on its own. 
In principle, it is easier to calculate the electrostatic self-energy of a point charge from the 
electrostatic potential in closed form rather than given in the form of series of multipoles. It 
may be practical to have an expression easy to handle for some calculations. 

This paper is devoted in the resolution of the electrostatic potential in a particularly simple 
wormhole. The spacetime representing the simplest kind of static, spherically symmetric 
wormhole has been popularized by Morris and Thorne in a pedagogical paper [8] but it had 
been already proposed by Ellis [9] and Bronikov [TO]. The metric describing this spacetime 
can be written 

ds 2 = -dt 2 + dl 2 + {I 2 + w 2 )(dtf 2 + sin 2 tidip 2 ) (1) 

in the coordinate system (t, I, </?), with — oo < I < oo, where it; is a constant. The space- 
time of this wormhole is regular everywhere. The throat is located at I = 0. It admits two 
distinct, asymptotical spacetimes which are obtained in the limit where I — > ±oo in metric 
these are two Minkowskian spacetimes. One can find an overview of wormholes in the 
book of Visser [TT] . 

We consider a charged particle at rest in this spacetime, assuming that the electromagnetic 
field is sufficiently weak so that its gravitational effect is negligible. In this case, the 
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electromagnetic field satisfies the Maxwell equations in the background metric (H|). The 
charge density p of a point charge e located at I = Iq, ■& = i)q and ip = po is given by 

P (i, 0, y>) = -^=<y(z - i )6($ - ^o)<% - ^o) (2) 

where 3 is the determinant of metric ([T]). 

The electrostatic potential V, defined from the electric field Fqi by F^i = —diV, obeys the 
electrostatic equation 
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dP P + w 2 dl {P + w 2 
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sin i? d"& \ d"9 J sin i? <9y? : 



-47rp. (3) 



To obtain a physical solution to equation ([3]) , we must add the asymptotic condition that the 
associated electric field Fqi vanishes at the two spatial infinities, i.e. as I — > ±00, because 
only the electric field has a direct physical meaning. 

A solution to equation ([3]) was given for the first time in the form of series of multipoles 
by Khusnutdinov and Bakhmatov pQ. Recently, Beskin et al |5] took up the issue of the 
multipole expansion. 

However, one of us in a previous preprint [12] has solved the electrostatic equation ([3]) in 
isotropic coordinates giving an explicit expression of the potential V in closed form. A part of 
this paper is devoted to clarify the relation between these two forms of the solution. Firstly, 
we find again the multipole expansion of V in a simple and synthetic manner. Secondly, we 
enlarge the result of the preprint [12J to a wormhole in n + 1 dimensions, n > 3, described by a 
direct generalisation of metric ([1]). To facilitate the reading of the paper, we give the complete 
proof in an appendix. In the case n = 3, we express the explicit expression of V in terms of 
the radial coordinate I. This allows us to compare easily with the multipole expansion. 

We then prove that the two determinations of V coincide. However, this above mentioned 
potential V is a particular solution to equation ([3|). In fact, we can add the monopole solution 
which satisfies the boundary condition about the electric field at I —> ±00. Assuming that 
the electrostatic potential must be symmetric in I and Iq, the general solution V\ to equation 
([3]), physically acceptable, depends on an arbitrary parameter A. 

We also consider the spacetime resulting from an identification according to the Z2 sym- 
metry, described by metric ([T]) with I > 0. We determine the electrostatic potential Vz 
generated by a point charge in this case. In this case it is unique. 

We find the explicit expression W of the electrostatic self-energy for the cases V, V\ and 
Vz by using the usual procedure. By deriving W, we then deduce the electrostatic self-force 
which has only a radial component. 

The paper is organized as follows. In section 2, we give the expressions of the radial 
functions of the multipole solutions. We express the electrostatic potential in the form of 
series of multipoles in section 3. We discuss in section 4 the general potential which is 
physically acceptable. We determine in section 5 the electrostatic potential in closed form 
and we prove that it coincides with the above multipole expansion. We obtain the electrostatic 
potential in the wormhole with the Z2 symmetry in section 6. We perform the calculation of 
the electrostatic self-energy and of the self- force in section 7 in all cases considered. We add 
some concluding remarks in section 8. 
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2 Radial functions of the multipole solutions 



To simplify the writing of the formulas, we introduce the dimensionless variables x and a 

I 



w 



la 

w 



and the potential without dimension A by putting 



w , 



A(x,#,<p) = -V(x,ti,tp) 



(4) 



(5) 



The electrostatic equation (jSJ) becomes 
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We will seek the multipole solutions to the corresponding homogeneous equation 
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(6) 



(7) 



They are obtained by separation of variables x, i? and In spherical symmetry, it is well 
known that this is performed by using the spherical harmonics Y™($, (p), n = 0, 1, . . . with 
\m\ < n. The multipole solutions have thus the form 



AZ{x,0,v)=A n {x)Y?(ti,ip). 



(8) 



By substituting form (|SJ) into equation (J7J), we find that each radial function A n satisfies the 
differential equation 

d 2 A n + 2x dA n _ n(n + l)A n = Q 
dx 2 x 2 + 1 cix x 2 + 1 

The coefficients of the linear differential equation Q are regular everywhere and therefore 
the vector space of the solutions consists of smooth functions. The solutions have two types of 
asymptotic behavior: x n and l/x n+1 as x — > oo. We have to define two linearly independent 
solutions corresponding to the two asymptotic behaviors that we denote respectively f n and 
g n , the latter been fixed up to a multiplicative constant. 

We introduce the variable z by setting 



IX. 



or x 



-iz. 



where % is the unit imaginary number. The differential equation ([9]) expressed in z becomes 
the Legendre differential equation whose the theory in the complex plane is well known |13j . 
The two fundamental solutions in the complex plane are denoted P n (z) and Q n {z). The P n (z) 
are polynomials in z of degree n whose powers of z are either odd or even. The functions 
Q n ( z ) can be written 

Quiz) = -P„0)hi +w n -i(z) 
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where the w n -i(z) are polynomials of degree n — 1, with W-\ = 0, whose powers of z are 
either odd or even, given for instance by 



„ 1 



Wn-l{z) = -Pq-l(z)P n -q(z). 
q=l 1 



The Wronskian of these solutions is 

W(P n ,Q r 



z 2 - 1' 

They admit respectively the asymptotic behaviors 

2 n T{n + \)z n 

p n{z) ' 7 == r , — —pr- as z -> oo, 

V7rr(n + 1) 

„ , , v^r(n + 1) 

Qn(z) ~ — n — ; 51 rr as z — )• oo. 

; 2 n + l T(n + f)^ 1 

The above results suggest immediately that one can take 

f n (x) = ^Pn(ix) (10) 

as a real solution to the differential equation Q with the asymptotic behavior 

, 2 n T{n+\)x n 

fn{x) J==- — pr- as X Y OO. (11) 

V7rr(n + 1) 

On the other hand, we cannot take Q n (ix)/i n as the other solution g n because the 
logarithmic term in Q n {ix) yields a problem at x = 0. Indeed, the expression of Q n (z) for 
z = ix is 

I fix + 1\ 

Qn(ix) = ~ P n(lx) In + W n -i(ix). 

I \ix — 1 J 

But we have the identity 

1 j, (%x + l\ 

arctan — = — In , x f 0, 

x 2 \ix — 1/ 

and therefore 

Q n (ix) = arctan — h w n _i(zx), if 0. 

z x 

There exists a discontinuity of the function Q n (ix) at x = 0. 

Considering the above expression, we define <7 n (x) for positive x by putting 

9n(x) = (-l) n 90-, x>0 . (12) 

The factor (— l) n is introduced for later so that the Wronskian of f n and g n does not depend 
on n. Taking into account the identity 

1 vr 

arctan — = arctan x, x > 0, 

x 2 
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we write g n (x) for x > under the form 



g n (x) = (-iy 



f n (x) - arctanx^ + v n -i(x) 



(13) 



where the v n -\(x) = w n -i(ix)/i n ~ 1 are real polynomials in x of degree n — 1. 

But of course, we can now extend expression (|13p for x < 0. Thus, expression Q13]) is a 
real solution to the differential equation whose asymptotic behavior is 

\/vrr(n + 1) 

<?n (a? ~ ; — ; 7 as x — > oo. (14 

^ > 2 n + 1 Y{n+l)x n + 1 V ' 

Taking into account the identity 

1 vr 

arctan — = arctanx, x < 0, 

x 2 ' 

we note that an equivalent way to express (|13p for negative x is the formula 

5n(x) = (-ir(^^+vr^), *<0. (15) 

We have therefore found two regular solutions f n and g n to the differential equation ([9]) 
which are linearly independent, with the Wronskian 

W(/ n , 5 „) = -— (16) 
1 + x z 

and which have respectively the asymptotic behaviors (jlip and ()14p as x — >• oo. 

We may also consider the asymptotic behavior as x — > — oo of the solutions to the dif- 
ferential equation ([9]). The form of this latter shows that starting from a solution we obtain 
a new solution by remplacing x by —x. So we can define the solution g n (—x). It has the 
asymptotic behavior l/(— x) n+l as x — > — oo. 

By using the parity properties of polynomials /„ and v n —i, it can be easily shown that 

g n (-x) = irf n (x)-(-l) n g n (x), (17) 

and according to (|16p we have 

W(gn(-x),g n (x)) = -—^. (18) 

1 + x z 



3 Electrostatic potential in series of multipoles 

The form of the multipole solutions (|HJ) suggests to seek the solution A to the electrostatic 
equation (|6|) in the form of series of multipoles. By expressing the second member of equation 
([6]) in terms of spherical harmonics by the well known formula 

oo n 

5(co S tf-co S? V>^-¥>o) = £ E Y ?*(0o,<Po)Y?(0,<p), (19) 

n=0 m=—n 
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we see easily that the potential A has the form 

oo n 

A{x,$,y) = Y J E A n (x,a)Y™*(ti , Vo )Y™(ti,<p), (20) 



n=0 m- 



where * denotes the complex conjugate. By reporting (|20p into ([6]), we find that each com- 
ponent A n has to satisfy the differential equation 

d 2 A n | 2x dA n n{n + l)A n _ 4tt ^ 
dx 2 x 2 + 1 dx x 2 + 1 1 + x 2 



To solve equation (I2ip we must add the asymptotic behaviors of A n as x — > ±oo. As 
stated in the introduction, we require that the asymptotic electric field vanishes. Taking into 
account the asymptotic behaviors of the radial functions, it is sufficient to require 

lim A n (x,a) = 0, lim A n (x, a) = 0, (22) 

X — >— oo x— >oo 

defining thus a solution A to the electrostatic equation ©. 

We now use the results about the Green function of a self-adjoint differential equation, 
subject to asymptotic behaviors (|22p . given in the appendix A. We apply for the case 

p(x) = 1 + x 2 , q(x) = —n(n +1), c = —Att. 

The homogeneous equation associated to equation (f2T|) is identical to equation ([9]) which ad- 
mits the two regular, linearly independent solutions g n {—x) and g n (x), verifying the asymp- 
totic behaviors 

lim g n (-x) = 0, lim g n (x) = (23) 

x — y — oo x—>oo 



as shown in the previous section. Taking into account (|T8|) , we put 

A®{x) = 2g n {-x), A ( t\x) = 2g n {x). (24) 

We thus rewrite the Green function HA.Th in terms of the functions g n in the form 

A n (x, a) = A6(a - x)g n (-x)g n (a) + 49(x - a)g n (-a)g n (x), (25) 

where 6 is the Heaviside step function. 

To see that our result coincides with the one of Khusnutdinov and Bakhmatov [I], we must 
replace our radial functions g n (x) and g n {—x) by their expressions in terms of the Legendre 
functions according to the sign of x. For g n {x), we have (fl~2|) and (fl~5]) and from ([T7J) we get 

gn(-x) = 7T— X > 0, (26) 

9n(-x) = -90l, x<0 . (27) 

It will be useful later to know the asymtotic behavior of the potential A as I — > oo, 
or equivalently the electric flux $(oo) & t the infinity / —> oo. According to (|14p and (|25|) . 
the leading term in 1/x in expression (|20h of the potential A can only come from the term 
4go( — a)go(x). Let us note that 

7T 

/o(x) = 1, 5o(aj) = 2 ~ ar ctanx, (28) 
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and that Yq = \j\f^n. So we have 

4 1 fir \ 
A(x,$, ip) ~ — i - 7: + arctan a) as x — > oo. (29) 
X 47T \ 2 / 

Expression (I29p yields the electric flux $(oo) & t the infinity / — > oo given by 

$ (oo ) = 47re - + - arctan — . (30) 
v \2 i it)/ 

4 Physically acceptable general solution 



Solution (l25j) allows us to write an electrostatic potential (1201) satisfying the boundary con- 
dition about the associated electric field. This is not the general solution to the electrostatic 
equation ([6]) because we can add some homogeneous solutions A™ given by ([8]). However, 
only the electric field associated to the monopole solution Aq, given by 

A° (x) = -t=/o(s), A° (x) = -j=9o(x), 



vanishes both at the two spatial infinities, i.e. x —> ±oo. 

In consequence, the physically acceptable general solution to the electrostatic equation 
([6]) has the form 

A gen (x, if) = A(x, ■!?, if) + X(a)g (x), (31) 

where A is an arbitrary function of a. It cannot be determined by the study of the electrostatic 
equation (|6|) without another condition. Since the source term is symmetric in x and a, we 
may take the assumption that the general solution must be symmetric in x and a. This 
determines the function A in the form A(a) = Xgo(a) where A is now an arbitrary constant. 
So we put 

A x (s, i?, ip) = A(x, ip) + Xg (a)g {x). (32) 

For the general solution Ax given by (j32]h we obtain the electric flux $a(oo) by adding the 
monopole contribution to the electric flux ([30]) . We find 

&x(oo) = ^ 7Te ( — 1 — arctan — ^ + 47reA ( arctan —\ . (33) 

y ' \2 TT w J \2 w J 

Beskin et al [5] suggested the need to impose the additional condition that the electric 
flux at the infinity I — > oo is exactly Aire as expected in the usual cases at the spatial infinity. 
From expression (|33p . we obtain this condition for A = \jix. However, this choice favors one 
of the two infinities. 



5 Electrostatic potential in closed form 

We give in appendix B the electrostatic potential in a wormhole described by metric (pQ) 
written in n + 1 dimensions, n > 3. Its explicit expression is found in isotropic coordinates 
of this metric. In this section, we recall only the basic formulas of appendix B useful to 
understand the expression of the electrostatic potential V. 
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We firstly write metric ([I]) in isotropic form by introducing the radial coordinate r by 

,2 



w- 



1 W 

or I = r , 

4r 



and then metric fl]) takes the form 



(34) 



ds 2 = -dt 2 + 1 + ^2 ) ( dr2 + r2 ( d $ 2 + sin2 W 



(35) 



with r > 0. The position ro of the point charge is given by 



or / = r - - — . 

4r 



(36) 



Due to the spherical symmetry, without lost of generality, we write down the expression 
of the electrostatic potential in metric (j35[) for a point charge e located at r = tq and i?o = 
by using (|B.8P with (|B.12p for n = 3. We have the form 



V(r, i?) 



F(s) 



r + 



4r 



w 



(37) 



where s defined by (|B.14|) reduces to 

w 2 (r 2 — 2rr$ cos $ + Tq) 



4(r 2 + u; 2 /4)(r 2 + w 2 /4) ' 
and where F(s) is given by (|B.2ip for n = 3 in the form 



4r / 



< s < 1, 



2^^(^ 



7T 



1 H — arcsin(l — 2s) 



7T 



(38) 



(39) 



Expression f)37j) with (J39J) yields an electrostatic potential V which is regular everywhere 
except at the position of the point charge, characterized by s = 0. 

In this paper, we add the determination of the electrostatic potential (|37h in terms of I 
and Iq. By substituting (|34"|) and ([36]) into (|37|) . we have firstly 



1 - 2s 



ZZn + W 2 COS $ 



\// 2 + w 2 Jl 2 l + w 2 



(40) 



and finally we obtain 



^l 2 - 2ll cos $ + ll + w 2 sin 2 



1 1 

— I arcsin 

2 7T 



VZ 2 + w 2 \An + ^ 



(41) 



Of course, potential (1411) is a solution to electrostatic equation ([3]). The associated electric 
field vanishes as I — > ±oo. We point out that expression (|4ip of the electrostatic potential V 
is much simpler in these coordinates. 
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We now determine the electric flux ^(oo) as / — )• oo corresponding to solution (|4ip . We 
find immediately 



*(«,) = 47re 



1 1 

— I arcsin 

2 7T 



With the help of the identity 



(42) 



arctan a = arcsin 



we see that the electric flux (I42p coincides with the one of the solution in series of multipoles 
given by (|30j) . 

Since the only possible difference between two physical solutions is proportional to the 
monopole solution go(l/w), we therefore conclude that the potential in closed form and the 
previous potential in series of multipoles can be identified. Of course, we may add to V the 
monopole solution for obtaining the general solution V\ which is physically acceptable 



V x (l,0,<p) = V(l,0,<p) + X-g (lo/w)g (l/w). 

w 



(43) 



6 Wormhole metric with the symmetry 

We have considered until now that the spacetime is described by metric (pQ) with — oo < / < oo. 
However, we can impose the Z2 symmetry by the topological identification of the points 
(1,$, <p) and (—l,"&,<p) of this spacetime. Such a Z2 symmetry for massive wormhole models 
have been proposed as black hole foils [H] . 

The metric describing this new spacetime is again metric (H|) but with I > 0. However, the 
physical laws in the background metric (pQ) must satisfy some boundary conditions at I = 
compatible with the Z2 symmetry. In electrostatics, as already noticed in |14j . the radial 
component of the electric field must vanish at I = 0. 

We take a point charge e located at I = Iq and i? = with Iq > 0. The electrostatic 
potential Vz satisfies the electrostatic equation ([3]) for I > with the boundary condition 



dV z 



01 



0. 



1=0 



(44) 



It is easy to determine Vz by introducing, before the topological identification, a point 
charge e located at I = — Iq and 1? = 0. We get from (|4T|) the expression 



Vz{l,#) 



I 2 - 2ll cos & + Iq + w 2 sin 2 ■d 



11 I ll + w 2 cos d 

— I arcsin — 

2 7T \ y J l 2 +w 2 h2 + w 2 



+ 



I 2 + 2/Zn cos <& + ll + w 2 sin 2 1? 



1 1 

— I arcsin 

2 7T 



-Uq + w 2 cos 7? 



Vl 2 + w 2 Jl 2 + w 2 



(45) 



which does satisfy condition (I44p . We emphasize that solution (|45l) is unique because we 
cannot add the monopole solution go(l/w) which does not verify the boundary condition 
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By taking the limit of expression (|45p as I — > oo, we can calculate the electric flux $z(oo) 
at infinity. We find thereby 

$z(oo) = 4vre, (46) 

as expected from the Gauss theorem since there is no electric flux through the sphere I = 
by virtue of 



7 Electrostatic self-force 

Besides a divergent electric field at the position of the point charge, there is a finite electric field 
which may be considered as an external field. It exerts on the point charge an electrostatic 
force. So, we have firstly to exhibit the singular part of the potential which leads to the 
divergent electric field. 

We start with the electrostatic potential (|37p written in isotropic coordinates. According 
to (IB.19j) for n = 3, the singular part V s i n g of this potential is simply 

e 1 

v/l + ^/^V 1 + w2 / 4r o [ (x ' 

Since the function F given by (|39|) admits the expansion 

„, „ ew ( 1 2 1 ~-\ . . 

F(s) ~ — [—= + -Vs as s -> 0, (48) 



2 Vv 7 ^ t 2 
and taking into account (|37p . we obtain 

V(a?) ~ V sing (x) t r^-^ ^- as x -)■ x . (49) 

The regular term in the neighbourhood of xq can be seen as resulting from the difference 
between V and the elementary solution in the Hadamard sense to the electrostatic equa- 
tion. The method of regularization from the electromagnetic Green function in an ultrastatic 
spacetime leads to simply subtract V s i ng to V [E] . 

So the electrostatic self-force derives from the self-energy W as follows 

dW 



dx l 



in which the self-energy W is given by the formula 



W(x ) = - (V{x) - V sing (x)) as x -»■ x . (50) 



According to (|50p . we find from formula ([49 P the explicit expression 



W(r ) = - tT\,AV ^ 
27r(rQ + w z /A) z 



The electrostatic self-energy (|5ip can be written with the radial coordinate I in the form 



9 

ew 
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Result (|52p coincides with the one of Khusnutdinov and Bakhmatov [T] found within the 
multipole formalism. 

We now consider the general electrostatic potential (|43|) . Since the monopole contribution 
is symmetric in I and Iq, we may use formula (|50p . By adding to self-energy (|52p the monopole 
contribution, we find the self-energy W\ 

, , e 2 w Xe 2 fir l \ 2 , . 

W\(Iq) = -—-2 ^ + — - - arctan - . 53 

In the spacetime resulting from an identification according to the Z<i symmetry, we obtain 
from (|45|) the expression of the self-energy Wz for Iq > 



e 2 -w e 2 



Wz{l0) = - 2*®+v?) + 4/, 



1 1 . f-l 2 + w r 
- + — arcsm —k — 

2 vr Wq + w 



(54) 



We point out that the self-energy Wz given by ([54"]) becomes large when Zo is close to zero. 

We now discuss the electrostatic self- force. It has only a radial component /. From 
self-energy ([52]) corresponding to the potential V, we get 

•f( /o ) = fit Wl °2\2 - ( 55 ) 

We see that self- force (|55p is always attractive with respect to the throat. It vanishes at I = 0. 
For the general potential V\, we obtain immediately from self-energy (|53j) 

f^( l o) = in Wl °2\2 + 2 f ~ arctan -) • ( 56 ) 

If we adopt the value A = 1/ir which ensures that the electric flux at infinity is 4-7re, we find 
that the radial component fi/ w is always positive. 

In the spacetime considered with Z<i symmetry, we obtain from the self-energy (|54p . for 
Zo > 0, 



,. , e 2 wl e 2 e 2 w e 2 . { -In + w 2 \ 

fM = -,(i 2 + w 2 ) 2 + W + + i^g arcsm l-fiTtf I • ( 57 ) 



The self-force fz is always repulsive with respect to the sphere I = 0, edge of the spacetime. 
It diverges with the leading term 1/Zq as Iq — > 0. 



8 Conclusion 

We have found the physically acceptable general solution V\ to the electrostatic equation ([3]) 
in closed form, given by (j4"3j) with (|4"T|) . It depends on a parameter A due to the fact that 
the electric field of the monopole solution is regular both at the two spatial infinities. This 
arbitrariness is physically unsatisfactory. We may choose A = 1/ir so that the electric flux 
^A(cx>) a t the infinity I — > oo is equal to 47re. Infinity at I — > oo is preferred in this case 
and we think that this choice is physically contrived. From this point of view, the wormhole 
considered with Z<i symmetry is interesting since there is unicity of the electrostatic potential 
Vz with obviously <&z(<x>) = 47re. 
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We have performed the calculation of the electrostatic self-force in all cases considered. 
For the potential V, the self-force is always attractive with respect to the throat I = 0. 
On the other hand for the potential V\ with A = 1/-7T, we have found a self- force always 
directed toward infinity I — > oo. This result is radically different from the previous one. For 
the potential Vz in the spacetime considered with Z2 symmetry, we have found a self-force 
always repulsive with respect to the throat, edge of the spacetime. It diverges when the point 
charge tends to the throat. Consequently, the charge can be held at rest in this case only 
with an external force increasingly large. 



A The Green function of the radial differential equations 

Let the self- adjoint differential equation with a source term 

^ p{x)^-\ + q(x)A = c5(x - a), (A.l) 



dx \ dx 
and the corresponding homogeneous equation 

£{p( x )i-) + q( x )A = 0, (A.2) 
dx \ dx J 

where the functions p, dp/dx and q are continuous and moreover p is strictly positive on 
(—00, 00). 

In the context of our present work, we assume that the homogeneous equation (|A.2p . 
admits two regular solutions A® and A^ e \ linearly independent, having respectively the 
following boundary conditions: 

lim A®(x) = 0, lim A^ix) = 0. (A.3) 

The solutions A^ and A^ are fixed up to a multiplicative constant. The nonvanishing 
Wronskian W{A^,A^) can be easily expressed in the form 

W(A^,A^)(x) = -^-, (A.4) 
p{x) 

where k is a constant. 

The Green function A(x,a) of the differential equation (jA.ip . which verifies the boundary 
conditions 

lim A(x,a) = 0, lim A(x, a) = 0, (A. 5) 

x— >— OO x— >oo 

is defined by the other following properties: 

a. A(x, a) satisfies the homogeneous equation (|A.2|) except at x = a. 

b. A(x, a) is continuous at x = a. 

c. The derivative of A(x, a) with respect to x has a jump discontinuity of magnitude c/p(a) 
at x = a. 
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These conditions ensure that A(x, a) satisfies the differential equation (|A.1|) in the sense of 
distributions with the boundary conditions (|A.5|) . 

It is now possible to determine this Green function A(x, a) with the help of the functions 
AW and A^ e \ Consider the function H defined by 

c f A^{x)A ( - e \a) if 
H[X ' a) = k\A^{a)A^{x) if !> fl * (A - 5) 

It is obvious that H (x, a) verifies all the properties of a Green function. For simplicity, 
it is possible to take k = c by a judicious choice of functions A® and A^ defined up to 
a multiplicative constant. Thus, the Green function of the differential equation (|A.1|) has 
expression (1A.6|) with k = c that one can rewrite in the form 

A(x, a) = 9{a - x) A®(x) A^ e \a) + 9(x - a)A^(a)A^(x), (A.7) 

where 9 is the Heaviside step function. 

B Explicit electrostatic potential in higher dimensions 

Metric ([I]) can be easily generalised to describe the spacetime representing a wormhole in 
n + 1 dimensions, n > 3. We take 

ds 2 = -dt 2 + di 2 + (i 2 + w^dn 2 ^, (b.i) 

with — oo < I < oo, and where dVt^^ is the metric of the sphere S n ~ l described in the 
coordinates . . . , i9 n _2, </?)• The change of radial coordinates (|3"1|) expresses metric (jB.ip 
in isotropic coordinates 

2 

/ ;/•"• \ 



* 2 +(l + ^) (dr'+r^nU) (B.2) 



with r > 0. We set the associated Cartesian coordinates (x l ), i = 1,. . . ,n. Metric (|B.2p can 
be written thereby 



ds 2 = -dt 2 + (l + ^\ ((dx 1 ) 2 + --- + {dx n f), (B.3) 



where r = yj (x 1 ) 2 + • • • + (x n ) 2 . 

We consider the usual Maxwell equations in a spacetime of n + 1 dimensions. The charge 
density p for a point charge e located at X — Xq has the expression 

p{x i ) = ^/=<5 (n) (x* - 4), (B.4) 
where g is the determinant of metric (|B.3[) given by 

/ 9 \ n 
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(B.8) 



The electrostatic equation in metric (|B.3|) takes the form 

r* W 2 e , i • 

+ fc(r) 7 fl,V = -2(n - 2) f ^y- ; ^=j* ( *V - 4), (B.6) 

where A is the usual Laplacian operator, and in which the function h has the expression 

1 d f w 2 Y l ~ 2 _ (n-2)w 2 

[r) / w 2\ n - 2 dr[ 4r 2 ) 2r(r 2 +w 2 /4)' ( ' 

y + ¥2 J 

In a previous work |16j , it was shown that exact solutions to the electrostatic equation 
(|B.6P exist in the form 

with 

r(x, x ) = (x 1 - xl) 2 + ... + {x n - x%) 2 , (B.9) 

if the function h satisfies a certain differential equation. It may be only the type II with 
a = w/2 in the classification of this work. Consequently, the function h must satisfy the 
differential equation 

d, 1. n-1 2[-n(n-2)w 2 /4- Ab 2 } 

-rh + -h 2 + h--± V . 9 ' „' - = B.10 

dr 2 r (r 2 + w 2 /4) 2 

in which the two constants b et A are to be fixed. 
From (lB~7j) . we calculate 

d i 1 , 2 n-l 1 (n 2 -6n + 8)w 2 

—h + -h 2 -\ h = -^— ,\ . (B.ll) 

dr 2 r 2(r 2 +w 2 /4) 2 y ' 

By comparing (IB.lOj) and (jB.llj) . we find the value of the constants: b = w/2 and A = — 4n+8. 

We are thus in the case where the electrostatic equation (|B,6P has solutions in the form 
(IB.8jl . According to the results of the work [16] , the functions k and g are expressed as 

r 2 + w 2 /4 r «-2 
* (r) = - S/T"' g < r)= (,^ + a ,V4)^ - <R12) 

and the function F is determined by the differential equation 

4(-s 2 + s)F" + (-4ns + n)F' + (-4n + 8)F = 0. (B.13) 

As variable s we have set ^ 

k(r)k(ro) ^ 
The value s = corresponds to the position of the point charge. We emphasize that 

< s < 1. (B.15) 
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This comes from that fact that we can always take the particular case where the point charge 
is located at 0i = without lost of generality, due to the spherical symmetry. Inequality 
()B.15p is then written 

w 2 (r 2 - 2rr cos-i?i + r 2 ) < 4(r 2 + w 2 /4)(r 2 l + w 2 /A), 

which is equivalent to the recognized inequality 

4r 2 rQ + 2rrQW 2 cos 01 + w 4 /A > 0. 

We firstly find a particular solution to (|B.13j) 

A classical method shows that the general solution to (1B.13|) is given by 

F(s) = Ri i fl + h Jo lu{1 ~ u)r '" du } (B ' 17) 

where f\ and f% are two arbitrary constants. The integral in (|B.17p can be explicitly performed 
for all values of n. 

We will choose the constants f\ and fi first of all so that solution (|B.8|) is regular every- 
where, except at x = xq. This occurs if the function F, given by (|B.17h . has only a singularity 
at s = 0. Now we have 

[u(l - u)] n l 2 - 2 du - f 1 [u(l - u)] n/2 - 2 du ~ (1 - s) n/2 - 1 as a -> 1, 

therefore we take the relation 

fi + h f 1 Hl-u)] n / 2 - 2 du = 0. (B.18) 
Jo 

The important point is that solution (|B.8|) satisfies the electrostatic equation (|B.6[) in the 
sense of distributions. By noticing from (1B.17|1 that 

we can deduce the most divergent part of V as x — > xq 

( 2\ n ~ 2 1 1 

V(x) ~ /! J- J - - ^ /4r 2 )n/2 -l (1 + w 2 / 4r 2)n/2-l jf^T x )]™/2-l ' 

By considering the source term of equation (|B.6p . we take 

/l = e (f) n_2 ' (R20) 

since 

1 ir n / 2 

A F(^F^ = - 2( "- 2) fw2) ,( " )fe ^»- 

So ,we obtain 

F(s) = e f-V 2 — * /0 = { 1 - ^ ; f'[u(l - u)] n l 2 - 2 du\ . (B.21) 

Now, with (|B.12|) . (|B.14p and (|B.2ip we are in position to write from (|B.8p .the potential V(x) 
generated by a point charge at x = xq. 
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